Orbital Polarons in the Metal-Insulator Transition of Manganites 
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The metal-insulator transition in manganites is strongly influenced by the concentration of holes 
present in the system. Based upon an orbitally degenerate Mott-Hubbard model we analyze two 
possible localization scenarios to account for this doping dependence: First, we rule out that the 
transition is initiated by a disorder-order crossover in the orbital sector, showing that its effect on 
charge mobility is only small. Second, we introduce the idea of orbital polarons originating from 
a strong polarization of orbitals in the vicinity of holes. Considering this direct coupling between 
charge and orbital degree of freedom in addition to lattice effects we are able to explain well the 
phase diagram of manganites for low and intermediate hole concentrations. 

PACS number(s): 72.80.Ga, 71.30.+h, 71.38.+i, 75.50.Cc 



The doping dependence of the properties of mangan- 
ites, .Ri-ajAcMnOs (where R and A represent rare-earth 
and divalent metal ions, respectively), poses some of the 
most interesting open problems in the physics of these 
compounds. First to be noticed is the peculiar asym- 
metry of the phase diagram that is most pronounced in 
the charge sector: Regions of high (x > 0.5) and low 
(x < 0.5) concentration of holes are characterized by such 
contrasting phenomena as charge ordering and metalic- 
ity, respectively. In the latter region - which we wish to 
focus on - the metallic state can be turned into an in- 
sulating one by raising the temperature above Tq » 200 
- 350 K. Introducing the notion of double exchange |l]] 
which associates the relative orientation of localized Mn 
tig spins with the mobility of itinerant e g electron, early 
work has identified this transition to be controlled by 
the loss of ferromagnetic order inherent to the metallic 
state. It is believed that the lattice effects are also of 
crucial importance in this transition. Within the lattice 
polaron - double exchange picture Q , the crossover from 
metallic to insulating behavior is controlled by the ratio 
of polaron binding energy to the kinetic energy of charge 
carriers. Assuming the former quantity to be constant, 
a critical coupling strength leading to localization can be 
reached by reducing the kinetic energy, i.e., by raising 
temperature which acts via the double-exchange mecha- 
nism. Spin disorder and spin polaron effects H further 
enhance the carrier localization above Tq- The doping 
dependence of the metal-insulator transition, however, 
is not fully captured in this picture. Namely the com- 
plete breakdown of metalicity at hole concentrations be- 
low x cr it ~ 0.15 - 0.2 that occurs despite the fact that the 
fully saturated ferromagnetic state is sustained remains 
an open problem which we will address in this paper. 

We analyze two mechanisms that could drive the lo- 
calization of doped holes at small x. First, we explore 
the possibility of the metal-insulator transition to be con- 
trolled by a disorder-order crossover in the e ff -orbital sec- 
tor. We outrule this scenario by showing that the mo- 
bility of doped holes depends only little on the degree of 



order in the orbital background. Second, we introduce 
the concept of orbital polarons in an orbitally degenerate 
Mott-Hubbard system - this constitutes the main idea 
of the paper. Similar to spin polarons in correlated spin 
systems such as cuprates, the orbital polaron is a natural 
consequence of strong correlations and the double degen- 
eracy of on-site levels, the latter being a consequence of 
the e g -orbital degeneracy present in ferromagnetic man- 
ganites. We argue that the coupling between charge car- 
riers and orbital degrees of freedom yields a strong po- 
larization of e g orbitals in the neighborhood of holes that 
acts together with breathing-mode like lattice distortions 
to form combined orbital-lattice polarons. The effect is 
more pronounced at small x because the polarizability 
of orbitals enhances as the characteristic energy scale of 
orbital fluctuations oc xt reduces. Within the orbital- 
lattice polaron scenario we are able to explain well the 
phase diagram of manganites for low and intermediate 
doping levels. 

Orbital order-disorder and hole mobility. — As the 
frustrating effect of holes is weakened at small x, a Jahn- 
Teller- and superexchange-driven transition from a dis- 
ordered, strongly fluctuating orbital state to an orbitally 
ordered state is expected. This possibly suggests the lo- 
calization process to be triggered by the establishment 
of order in the orbital sector, which would naturally in- 
troduce a doping dependence due to the energy scale xt 
of orbital fluctuations. We analyze this possibility by 
comparing the mobility of holes moving in orbitally dis- 
ordered and ordered backgrounds. The relevant transfer 
Hamiltonian of e g electrons has to account for strong on- 
site Coulomb repulsions and the orbital degeneracy: 
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with 7 6 {x,y, z}. We use constrained operators c\ a = 
cj (1 — m) that create electrons at site i in orbital a only 
under the condition that the site is empty. In orbital 
basis a € {\3z 2 — r 2 ),\x 2 — y 2 )} the intersite transfer 
matrices are: 
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Electron spins are treated within the double-exchange 
model and implicitly enter through the transfer ampli- 
tude t only. To observe the strongly correlated nature 
of Hamiltonian (|l|) it is convenient to introduce separate 
quasiparticles for charge and orbital degrees of freedom 
j|. In the case of an orbitally disordered state we em- 
ploy a slave-boson representation in which electrons are 
replaced by fermionic orbitons fj a carrying orbital pseu- 
dospin and bosonic holons b\ carrying charge 0]. On a 
mean- field level with parameters J2 a /3 t™ 13 (fiafjp)i = X 
and (b\) = ^/x these quasiparticles can be decoupled. 
The characteristic energy scale of orbital fluctuations is 
then set by the orbiton half band width -D orD = 3xt, the 
mobility of holes by the holon half band width Dhole = 3t. 

Orbital order in manganites is induced by a coupling 
of orbital pseudospins mediated by, e.g., J aim- Teller 
phonons and superexchange. The corresponding Hamil- 
tonian (see, e.g., is of XY-tyj>e with internal frus- 
tration making it difficult to handle. To simplify the dis- 
cussion we simulate an orbitally ordered state by adding 
to Hamiltonian (Q) the interaction term 
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6 and Pauli matrices erf. For q = (tv,tv,0) 

y 2 ))/V2 
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Eq. (||) favors a staggered (|3a 
ordering within x-y planes repeating itself along the z 
direction; this closely resembles the type of ordering ob- 
served experimentally in LaMn03 |tJ . The establishment 
of orbital order is indicated by a singularity in the static 
susceptibility (a x a x ) q . In random-phase approximation 
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with J q = J {cosq x + cosq y — cosq z ) /3. Bare suscep- 
tibilities (• • -)° are treated employing the above mean- 
field description. Numerically solving for the pole of Eq. 
(0) with J = 0.13 eV as estimated from the structural 
phase transition observed in LaMrrOs at T = 780 K |?]] 
and t = 0.4 eV we find a critical doping concentration 
£crit = J/4t = 0.08; at concentrations below this critical 
value an orbitally ordered state is to be expected. 

We now turn to analyze the effect of orbital order onto 
the mobility of holes. For this we study the evolution 
of the holon half band width, which in the disordered 
state -Dhoic = 3t, as orbital order develops. Foremost, an 
important difference between models with orbital pseu- 
dospin and conventional spin is to be noticed: In the 
latter systems spin is conserved when electrons are trans- 
ferred between sites. In contrast, the transfer Hamilto- 
nian (Fy) of the orbital model is nondiagonal in orbital 



pseudospin - the transfer matrices t"@ , ty@ , and can- 
not be diagonalized simultaneously. Hence there is a fi- 
nite amplitude for holes to hop without perturbing the 
orbitally ordered background, leading to the formation 
of a wide coherent hole band; for the specific type of or- 
der introduced above its half width is D^ e — 2t. This 
number seems to indicate a profound reduction in the 
mobility of charge carriers by ~ 30% as compared to the 
disordered state. It does, however, not reflect the full 
kinetic energy of holes which is further contributed to 
by incoherent processes involving the creation of excita- 
tions in the orbital background. To study the influence of 
these processes we employ an "orbital-wave" approxima- 
tion that uses a bosonic description of orbital excitations 
of energy J and a fermionic description of holons moving 
within a band of half width 2t. Holons are coupled to or- 
bital excitations via part of Hamiltonian (|l]) . In analogy 
to studies of spin systems || we analyze the motion of a 
single hole employing a self-consistent Born approxima- 
tion. Two limits can be treated analytically: As orbital 
excitation become very heavy, J/t — > oo, the hole decou- 
ples from the orbital subsector. It then moves coherently 
due to the nondiagonality of transfer matrices alone, its 
motion being maximally constrained as -Dhole = 2t. In 
the opposite limit, J/t — > 0, the hole moves completely 
incoherently in a cloud of an infinite number of orbital 
excitations suppressing orbital correlations in the back- 
ground. The hole mobility then reaches its maximum 
-Dhoic = 3t as in the disordered state, while the effective 
mass of the composite object formed of hole and excita- 
tion cloud becomes large, m — * oo. For a ratio J/t = 0.3 
realistic for manganites one is closer to the latter limit: 
Numerically we find the establishment of orbital order to 
reduce the kinetic energy of the hole by only less than 5%. 
We therefore conclude that a disorder-order crossover in 
the orbital sector has only a secondary effect on the mo- 
bility of charge carriers, outruling it as a driving mecha- 
nism to initiate the metal-insulator transition. 

Orbital polarons . — In the preceding paragraph we 
have considered charge carriers and orbitals to interact 
via the transfer part (|l]) of the Hamiltonian. We will 
now show that in an orbitally degenerate Mott-Hubbard 
system there also exists a direct coupling stemming from 
a strong polarization of e g orbitals on sites next to a 
hole. The lifting of orbital degeneracy is mediated by a 
displacement of oxygen atoms situated between the Mn 
sites as well as by the e g level splitting in the vicinity of 
positively charged holes. We estimate the magnitude of 
the degeneracy lifting, A = A ph + A ch , as follows: Holes 
respectively electrons are coupled to the lattice breathing 
mode Qi and to two Jahn- Teller modes Q2 and Q3 by 

#oi-ph = X) ( ~ 9lQi,i n i + 92 (Qa.iCTf + Qs.-tof) 
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where is the number operators for holes. Integrat- 
ing over oxygen displacements Q = {Qi,Q2,Q?,) we 
obtain A? h = gi g 2 V^/(3K) « (gi/g 2 )Ej T . A lower 
bound for this quantity is given by the Jahn- Teller en- 
ergy i?jT ~ 0.1 - 0.2 eV, assuming that coupling to the 
breathing mode is at least as strong as to the Jahn- Teller 
one. On the other hand, the contribution to the e g level 
splitting from the Coulomb interaction between the posi- 
tively charged hole and an electron on a neighboring site 
is estimated taking into account the covalency of Mn 3c? 
and O 2p orbitals. It follows A ch w f 7 2 e 2 i?Mn-Mn, where 
the covalency factor 7 = t p( i/ A p d can be obtained from 
the transfer amplitude and the charge gap between Mn 
and O sites, t p d ~ 1.8 eV and A p( j w 4.5 eV j|, respec- 
tively. With lattice spacing i?Mn-Mn = 3.9 A this leads 
to A ch ps 0.4 eV. In total, the polarization of e g levels 
on sites next to a hole yields an energy splitting A w 0.5 
- 0.6 eV. Being comparable in magnitude to the transfer 
amplitude t this number strongly indicates a direct cou- 
pling of charge and orbital degrees of freedom to be of 
importance in manganites. 

In the cubic system the Hamiltonian describing the 
above coupling is 
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where tJ = 



(sin Qa x + cos6cr z ) /2 denotes orbital pseu- 
dospin operators with angles = =f27t/3,0 correspond- 
ing to bonds in directions 7 = x,y, z. Hamiltonian (|5|) 
promotes the formation of orbital polarons. For low 
enough hole concentrations these consist of a bound state 
between a central hole and the surrounding e g orbitals 
pointing towards the hole. Besides minimizing interac- 
tion energy this configuration yields a large amplitude 
of virtual excursions of e g electrons onto the empty site, 
thus lowering the kinetic energy as well. In other words, 
the double exchange process is locally enhanced, provid- 
ing a large effective spin of the orbital polaron. This nat- 
urally explains the development of ferromagnetic clusters 
experimentally observed at temperature above Tq p0| - 
At large doping concentrations, orbital polarons cannot 
be considered as isolated anymore. In the presence of 
many holes, the fast fluctuations of orbitals destroy the 
orbital-hole bound state. To elaborate on this point we 
calculate the orbital-hole binding energy -D rb by consid- 
ering a static hole in an orbital background. The upper 
bound for -E or b is set by the limiting case of uncorrelated, 
static orbitals: Hamiltonian (j|) yields _E or b = 3A(1 — x) 
if all six orbitals surrounding a hole point towards the 
latter, where the factor (1 — x) accounts for the fact that 
some of the orbitals might be empty. In the real system, 
however, the binding energy is reduced by inter-orbital 
coupling with energy scale cx J and orbital fluctuations 
of frequency cx xt which are comprised in the effective 
"stiffness" of the orbital sector -D rb = 3xt + J. Calculat- 
ing the binding energy for a static hole in both orbitally 



ordered and disordered states we find an interpolating 
formula 
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which correctly reproduces the limits A 3> -D orD an d 
A -C -Dorb- In marked contrast to conventional lattice 
polaron theory, the binding energy of the orbital po- 
laron explicitly depends on the concentration of holes, 
foremost via the energy scale of orbital fluctuations xt: 
The more holes are present in the system the faster or- 
bitals fluctuate and the smaller the binding energy E m \, 
becomes. The crossover from the small-polaron regime 
in which the orbital-hole bound state is stable to the 
large-polaron regime in which the bound state breaks 
up is governed by the dimensionless coupling constant 
A rb = ^orb /-Dhole; where -D^oie is the half band width 
of holes. Since the lattice breathing mode in Eq. (Q) fur- 
ther enhances the polaron formation, the total coupling 
constant A = A or b + A p h has an additional contribution 
A p h = -B p h/-Dhoic, where E ph = gl/(2K). It is impor- 
tant to note that a critical value of A leading to local- 
ization can be reached either by increasing temperature 
which quenches -Dhole via double exchange or by lowering 
x which enhances -Sorb- 
in the following we study in more detail the conse- 
quences of the polarization of orbitals around doped holes 
for the phase diagram of manganites. We start from the 
strong-coupling limit A > 1 in which holes are substi- 
tuted for by small polarons that move in a band whose 
width is reduced by an exponential factor 
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where loq is the frequency of dispersionless phonons. 
In deriving Eq. (Q) we neglect the interference be- 
tween orbital and lattice couplings. To simulate the 
crossover from the small polaron to the free carrier pic- 
ture we phenomenologically employ the function 7 — [3 — 
In [A (1 + 13)] /A 2 with 13= (l- 1/A 2 ) 1/2 which has been 
proposed for strongly coupled phonon systems ||TT[| ; for 
A < 1 holes are assumed to move freely as 7 = 0. A feed- 
back of Eq. (0) into the binding energy, Eq. @, result- 
ing from the quenching of orbital fluctuations cx xte -71 
as small polarons begin to form sharpens the transition. 
Next we turn to the coupling between the charge and spin 
sector through the double-exchange mechanism. On the 

1 /2 

one hand, the transfer amplitude t = to [(l + m 2 ) /2] 
depends on the magnetization m. This is determined 
from m = Bs(am) with a = 35/ (1 + S)Tq/T, where Bs 
is the Brillouin function and S = | + (1 — x)h the aver- 
age spin of Mn 3c? electrons per site. On the other hand, 
charge transfer mediates a ferromagnetic interaction be- 
tween sites. The effective exchange coupling constant is 
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FIG. 1. Magnetic and electronic phase diagram obtained 
within the present theory. Dashed and dashed-dotted lines 
represent magnetic and electronic transitions, respectively, 
the simultaneous transition in both channels is denoted by 
a solid line. The phases are: paramagnetic insulator (PI), 
paramagnetic metal (PM), ferromagnetic metal (FM), and 
ferromagnetic insulator (FI). 



ferromagnetic metal to a paramagnetic insulator occurs. 
We note that the loss of charge mobility triggers static 
orbital order. A long-range orbitally ordered state has 
in fact been experimentally detected in the insulating re- 
gions of Lao.88Sro.i2Mn03 JlS} ]. In general, however, the 
ordered state is expected to have orbital and Jahn- Teller 
glass features due to the presence of quenched orbital po- 
larons. The orbital polaron provides a strong local ferro- 
magnetic coupling between the planes in c-direction, thus 
explaining the existence of the isotropic ferromagnetism 
even in the insulating phase at small doping. Finally it 
is worth to notice that the phase diagram in this theory 
is highly sensitive to the transfer amplitude to as this 
parameter enters in the polaron binding energy, Eq. (@). 

We conclude that orbital polarons are an intrinsic 
part of an orbitally degenerate Mott-Hubbard system 
and play an important role in the physics of mangan- 
ites. Being comparable in magnitude to lattice effects 
this new mechanism naturally introduces a doping de- 
pendence into the polaron binding energy via the degree 
of orbital fluctuations. Accounting for both orbital and 
lattice effects we are able to reproduce well the important 
features of the phase diagram of manganites. 



J cS = ^ [x(l - x)e-v + x (l - x) 2 ] (8) 

with xq = 2t\/U , where U is the on-site repulsion, and 
mean- field parameter \ ~ \- The Curie temperature 
is obtained from Tq = ^vzS{S + 1) J e g; the fitting pa- 
rameter v compensates for the overestimation of Tc in 
mean-field treatment. The first term in squared brackets 
of Eq. (||) stems from the coherent motion of polarons 
and represents the conventional double-exchange contri- 
bution to Tc- The second term is due to superexchange. 
It is insensitive to the localization process and dominates 
ferromagnetic interaction at low doping levels |l^ Q . 

To illustrate the interplay between the system of equa- 
tions (^) - (^) we numerically extract from them the T-x 
phase diagram. The following parameters are chosen: 
The orbital polarization energy is set to A = 0.5 eV, 
yielding a binding energy comparable to the phononic 
one -Eph = 0.6 eV; the phonon frequency is ujq = 0.05 eV, 
the interaction between orbitals J = 0.13 eV, the bare 
transfer amplitude to = 0.4 eV, and xq = 0.1. The fit- 
ting parameter v = 0.5 JlJ] is adjusted to reproduce the 
values of Tc observed for Lai^^Sr^MnOs [|l6|. The re- 
sult is shown in Fig. [j]. Due to the doping dependence 
of the orbital polaron binding energy the system is more 
insulating at low and more metallic at high doping levels. 
This is manifested in the absence of a metallic phase at 
doping levels x < 0.15 and the appearance of a paramag- 
netic metallic phase at x > 0.4. In the region in between 
a simultaneous magnetic and electronic transition from a 
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